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An overview is given on results from direct and indirect measurements of galactic cosmic 
rays. Their implications on the contemporary understanding of the origin of cosmic rays 
and the knee in their energy spectrum are discussed. 
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1. INTRODUCTION 

Cosmic rays (CRs) in the energy range from several GeV up to about 100 PeV are 
assumed to be mostly of galactic origin. At energies up to several 100 MeV indi- 
vidual isotopes can be identified, e.g. with the ACE/CRIS experiment, a satellite 
borne silicon detector telescope. At higher energies, CRs are identified by their 
charge and the energy measurement becomes an experimental challenge. Various 
techniques are utilized, like the determination of the particles momenta in magnet 
spectrometers (e.g. BES9=f), the (partial) absorption of nuclei in calorimeters (e.g. 



ATIC 3 ), or the measurement of transition radiation emitted by relativistic particles 
(e.g. TRACER^!) . Circumpolar long duration balloon flights offer the possibility of 
a long exposure (> 14 d) combined with low atmospheric overburden (typically 
< 5 g/cm 2 ) as recently demonstrated by the ATIG 5 , TIGER© and TRACEeP 
experiments. 

At energies above 1 PeV the steeply falling spectrum requires large detection ar- 
eas (exceeding several 10 4 m 2 ) and exposure times of several years, which presently 
can be realized only in ground based installations. They measure the secondary 
products generated by the CR particles in the atmosphere - the extensive air show- 
ers. The challenge of these investigations is to reveal the properties of the primary 
particle behind an absorber - the atmosphere - with a total thickness, correspond- 
ing to 11 hadronic interaction lengths or 30 radiation lengths. Consequently, these 
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experiments have a coarser resolution and only mass groups or the average primary 
mass are derived. Two basic approaches can be distinguished: Measuring the debris 
of the particle cascade at ground level by registering the main shower components, 
the electromagnetic, muonic, and hadronic parts. Or measuring the longitudinal 
shower development in the atmosphere, exploring the Cerenkov or fluorescence light 
generated predominantly by the shower electrons. Examples are the KASCADE - 
Grande^ or EAS-TOP^I installations, measuring simultaneously the electromag- 
netic, muonic, and hadronic shower components, the SPASE HI/ AM AND^Pl exper- 
iment, investigating electrons and high energy muons, or the BLANC^Pl Cerenkov 
detector. 

In this article an overview on results from direct and indirect measurements is 
given, concerning the sources of CRs their propagation through the Galaxy 
(©, and the energy spectra and mass composition observed at Earth (Q. Their 
implications on the understanding of the origin of the knee are discussed (jJHJ. 



2. SOURCES OF COSMIC RAYS 

A big step towards the understanding of CR sources would be their direct observa- 
tion in the sky. However, charged CRs are deflected in the galactic magnetic fields, 
the gyromagnetic radius of a proton with an energy of 1 PeV is about 0.4 pc. But 
7-rays, are good candidates for a point source search. Photon emission of supernova 
remnants (SNRs) has been detected in a wide energy range from radio wave lengths 
to x-rays. The observations are interpreted as synchrotron emission from electrons, 
which are accelerated in these regionJEl. The HEGRA experimental has detected 
an excess of 7-rays with TeV energies from the SNR Cassiopeia A. This is inter- 
preted as evidence for hadron acceleration in the SNR. The hadrons interact with 
protons of the interstellar medium close to the source region, producing 7r°s, which 
decay into high-energy photons. The flux is comp atib le with a model of electron 
and hadron acceleration in shock fronts of the SNEP-M 

Despite of the above-mentioned deflection, it is of great interest to study the 
arrival direction of ch arge d CRs as well. The result of such an analysis from the 
KASCADE experiment"*! is depicted in Fig. (left) . Shown is the distribution of 
the significances from a sky map of the arrival direction of showers with energies 
above 0.3 PeV covering a region from 10° to 80° in declination. For an isotropic 
distribution the significances are expected to follow a Gaussian distribution as indi- 
cated by the solid line. Results for all events are presented, as well as for a selection 
of muon-poor showers. The latter are expected from potential primary 7-rays. No 
significant deviation of the data from the Gaussian distribution can be recognized. 
The analysis has been deepened by investigating a narrow band (±1.5°) around the 
galactic plane. Also circular regions around SNRs and TeV-7-ray sources have been 
studied. None of the searches provided a hint for a point source. In addition, no 
clustering of the arrival direction for showers with primary energies above 80 PeV 
is visible. Claims by the MAKET-ANI experiment for a point-source detection^! 
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Fig. 1. Left: Distribution of the significance values from a sky map of the arrival direction of CRs 
as measured by the KASCADE experiments for the complete data set (open circles) and a selec- 
tion of muon poor showers (filled squares). Right: Comparison of derived CR source abundance s of 
refractory nuclides with solar-system abundances according to measurements with ACE/CRl£Hl 
normalized to 28 Si. 



have been withdrawn meanwhile"™ 

Despite no sources have been detected with charged particles, information on the 
composition at the source can be obtained from measurements of the abundance of 
refractory nuclei. They appear to have undergone minimal elemental fractionation 
relative to one another. The derived abundance a t th e source is presented in Fig.^ 
(right) versus the abundance in the solar systern^l The two samples exhibit an 
extreme similarity over a wide range. Of the 18 nuclides included in this comparison, 
only 58 Fe is found to have an abundance relative to 28 Si that differs by more than a 
factor of 1.5 from the solar-system value. When uncertainties are taken into account, 
all of the other abundances are consistent with being within 20% of the solar values. 
This indicates that CRs arc accelerated out of a sample of well mixed interstellar 
matter. 

Motivated by the observations, it is assumed that at least a large fraction of 
CRs are accelerated in supernova 

remnantPE™ 

However, recent progress in 

the understanding of 7-ray bursts has put forwar d the i dea that a subsample of 
high-energy CRs may be accelerated in 7-ray burst3^2E5J 



3. PROPAGATION OF COSMIC RAYS 

After acceleration, the particles propagate in a diffusive process through the Galaxy, 
being deflected many times by the randomly oriented magnetic fields (B ~ 3 fJ-G). 
The nuclei are not confined to the galactic disc, they propagate in the galactic halo 
as well. The diffuse 7-ray background, extending well above the disc, detected by the 
EGRET experiment, exhibits a structure in the GeV region, which is interpr eted 
as indication for the interaction of propag ating CRs with interstellar matteM 
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Fig. 2. Left: Me asured boron-to-carbon ratio as function of energy, the lines indicate model 
predictions, se Right: Rayleigh amplitudes as function of energy for v arious experiments, for 
references se f23 Additionally, model predictions for Leaky Box model and a diffusion modeC2 
are shown. The lines indicate the expected anisotropy for primary protons, iron nuclei, and all 
particles. 



The height of the halo has been estimated with measurements of the 10 Be/ 9 Be- 
ratio by the ISOMAX experimenJ23 to be a few kpc. The measured abundance of 
radioactive nuclei in CRs with the CMS instrument yields a residence time in the 
Galaxy of about 15 • 10 6 a for particles with GeV energieJZH. 

Information on the propagation pathlength of CRs is often derived from the 
measurement of the ratio of primary to secondary nuclei. The latter are produced 
through spallation during propagation in the Galaxy. As an exampl e, t he measured 
boron-to-carbon ratio is shown in Fig. [21 (left) as function of energy**! The energy 
dependence of the measured ratio is frequently explained in Leaky Box models by a 
decrease of the pathlength of CRs in the Galaxy A(R) = Aq(R/ Ro)~ 5 , with typical 
values A w 10 - 15 g/cm 2 , S w 0.5 - 0.6, and the rigidity R « 4 GV. 

At higher energies such measurements are not feasible due to the limited mass 
resolution of air shower experiments. However, at these energies the large scale 
anisotropy is expected to reveal properties of the CR propagation. The Rayleigh 
formalism is applied to the right ascension distribution of extensive air showers 
measured by KASCADePOI. No hints of anisotropy are visible in the right ascension 
distributions in the energy range from 0.7 to 6 PeV. This accounts for all showers, 
as well as for subsets containing showers induced by predominantly light or heavy 
primary nuclei. Upper limits are shown together with results from other experiments 
in Fig. [2 (right). It presents the Rayleigh amplitude as function of energy. The 
experimental results are compared to the anisotropy expected from calculations of 
the propagation of CRs in the Galaxy. The data reflect a trend predicted by a 
diffusion modef^l This indicates that leakage from the Galaxy and consequently 
a decreasing pathlength A(E) plays an important part during CR propagation at 
high energies and most likely, also for the origin of the knee. 

Leaky Box models are successful at GeV energies as discussed above. In the PeV 
regime, however, they seem to be faced with some difficulties. Two versions of a 
Leaky Box modeP^, with and without reacceleration, seem to be ruled out by the 
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Fig. 3. Energy spectra for elemental groups a) 
protons, b) helium, and c) iron. Open sym- 
bols give res ults of direct measurements, for 
references seeSsGEl Filled symbols represent 
data from air shower measurements: KASCADE 
electrons/muons interpreted with two interac- 
tion modelfill (preliminary) , KASCADE si ngle 
hadrons 3 ®, and EAS-Top electrons/muons 3 -*^. 
The data are com pared to calculations by 
Kalmykov et ajEHI (. . SveshnikovJSSl (___), 
and the Poly-Gonato mode. 36 ( — ). 



anisotropy measurements, see Fig- EI This relates to the extremely steep decrease of 
the pathlength A oc E~ - 6 , yielding a t Pe V energies unrealisticly small values for A. 
Even for a residual pathlength modePS at 1 PeV the pathlength would be smaller 
than the matter traversed along a straight line from the center of the Galaxy to the 
solar systerrPSl 



4. ENERGY SPECTRA AND MASS COMPOSITION 

At energies below 100 TeV the energy spectra of individual elements have been 
measured with detectors above the atmosphere. Examples for protons, helium and 
iron nuclei are compiled in Fig. EI The measured spectra can be described by pow- 
erlaws. For the iron spectrum at low energies the modulation due to the magnetic 
fields of the heliosphere causes the flux suppression. Actual experiments, like ATld^l 
and TRACER 7 , as well as the proposed ACCES9^2 space project are expected to 
improve the experimental situation in the region around 0.1 to 1 PeV, where large 
uncertainties are visible in the figure. More precise fluxes in this region would be 
valuable to intercalibrate air shower measurements. 

The elemental abundance at 1 TeV is presented in Fig. 0] (left) as function 
of nuclear charge number for elements up to nickel. The experimental status for 
the heavier elements is summarized in Fig. 0] (right). All stable elements of the 
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Nuclear chargeZ Nuclear charge Z 

Fig. 4. Left: Abundance of elements (Z < 28) in CrJSEID at 1 TeV. Right: Relative abun- 
dance of CR elements (Z > 28) normalized to Fc= 1 from various expe rim ents around 1 GeV/n. 
For references seJSHEl For comparison, abundances in the solar systerrlSl are presented as well, 
normalized to Si (Jeft) and to Fe (right). 

periodic table have been registered in CRs. In both panels the CR abundance is 
compared to the abundance in the solar systeirP=l normalized to silicon and iron, 
respectively. The overall similarity of the two samples of matter, already seen in 
Fig. 2] is reflected here on a coarser scale. 

At higher energies, many air shower experiments have reported fluxes for all 
particles. A compilation is presented in Fig. [3] (left). The energy scale of the in- 
dividual experiments has been slightly normalized (±10%) in order to match the 
flux with that obtained by direct measurement J^l A good agreement between the 
experiments in the reconstructed shape of the spectrum is evident. The knee at 
~ 4.5 PeV and a smaller structure at ~ 400 PeV, the second knee, are visible. 

Most valuable to reveal the origin of the knee are measurements of the energy 
spectra for individual elements or at least elemental groups. KASCADE studied 
the influence of different hadronic interaction models used in the simulations to 
interpret the datcP^I. Two sets of spectra, derived from the observation of the elec- 
tromagnetic and muonic air shower components, applying an unfolding procedure 
based on the Gold algorithm and using CORSIKA**" with the hadronic interaction 
models QGSJET and SIBYLL are compiled in Fig. [31 for three elemental groups. 
As can be seen in the figure, the fluxes depend on the model used. KASCADE 
emphasizes that, at present, there are systematic differences between measured and 
simulated observables which cause the ambiguities of the spectra. These conclusions 
apply in a similar way also to other experiments. A correct deconvolution of energy 
spectra requires a more precise knowledge of the hadronic interactions. 

Fig. |3 also shows the spectrum of primary protons, which ha s be en derived 
from the flux of unaccompanied hadrons measured by KASCADePBI. The spec- 
trum is compatible with the proton flux as obtained from the unfolding procedure 
when using the QGSJET model. The EAS-TOP experiment published two sets of 
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Fig. 5. Left: Normalized all-particle energy spectra from different experiments. The lines indicate 
the average all particle spectrum and the contribution of galactic CRs. The knee at ~ 4.5 PeV 
and the second knee at ~ 400 PeV~ 92 ■ E^ are indicated. Right: The average flux of the mea- 
surements (left) is represented by the data points. Additionally, spectra for elemental groups with 
the indicated charge number range according to a parameterization of the measurements are de- 
picted, including a proposed contribution of ultra- heavy elements (Z > 28), extrapolated from 
measurements at GeV energies ("?"). For details and references sea 



spectra with different assumptions about the contribution of protons and helium 
nuclei deri ved from the measurements of the electromagnetic and muonic shower 
components^!. The resulting fluxes are indicated by two squares per primary en- 
ergy. To guide the eye, the solid lines indicate power law spectra with a cut-off at 
Z ■ 4.5 PeV. 

The dash ed li nes represent calculations of energy spectra for nuclei accelerated 
in supernovae^l It is assumed that the particles are accelerated in a variety of 
supernovae populations, each having an individual maximum energy that can be 
attained during acceleration, which results in the bumpy structure of the obtained 
spectra. The dotted lines reflect calculations of the diffusive propagation of particles 
through the GalaxySS. The leakage of particles yields a rigidity dependent cut-off. 
Comparison with the data may suggest a qualitative understanding of the energy 
spectra. However, for a precise quantitative understanding, detailed investigations 
of the systematic errors of the measurements are necessary and the description of 
the interaction processes in the atmosphere needs to be improved. 

While the elemental abundance is relatively well known at low energies from di- 
rect measurements (see Fig.0J, at higher energies, air-shower experiments provide 
information on mass groups or on the average mass. Frequently, the mean logarith- 
mic mass (In A), defined as (In A) = ^rjlnAj, where r, is the relative fraction of 
nuclei with atomic mass number Ai, is used to characterize the composition. (In A) 
is often derived from the ratio of particles measured at ground level. For a primary 
proton more electrons and hadrons and fewer muons are registered as compared to 
an iron induced shower with the same energy. The data from many experiments 
are compiled in Fig. [fJl (left). They exhibit an increase of (In A) as function of en- 
ergy in the knee region. The increase is compatible with expectations, assuming a 
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Fig. 6. Mean logarithmic mass of CRs reconstructed from a) experiments measuring electrons, 
muons, and hadrons at ground leveP^ ^ 47 48 an< j j,) observations of the shower maximum 
interpreted with a modification of The lines indicate expectations according to the 

Poly Gonato modeEHJ 



cut-off behavior of the flux of individual elements as indicated in Fig. [3] by the solid 
line. The second class of experiments reconstructs the average depth of the shower 
maximum X max from the observation of Cerenkov and fluorescence light. Using 
the model QGSJET to derive the mean logarithmic mass from the data results 
in a light ma ss composition at high energies in contradiction to the findings just 
mentioned^. Introducing modifications to QGSJET, namely lowering the inelastic 
cross sections and slightly increasing the elasticity of hadronic interactions, this 
discrepancy can be reduceoP^ and the mean logarithmic mass rises as function of 
energy, see Fig. (right) . 

The average experimental values from both classes of air shower measurements 
presented in Fig. are shown as light grey area in Fig. It represents the mean 
value ±1 standard deviation. The dark grey area represents the results of direct 
measurements above the atmosphere. This experimental situation will be compared 
to predictions of various models in the next section. 

A different interpretation of the experimental results is given in Fig. (right) . 
The average of the flux values shown in the left panel is displayed by the data points. 
The spectra for e lem ental groups are presented according to a parameterization of 
the measurements""-! which corresponds to the solid lines in Figs. Eland® where the 
agreement with the data has been discussed. Also shown is a proposed contribution 
of ultra- heavy elements (Z > 28), extrapolated from measurements at GeV energies. 
The individual spectra exhibit a cut-off at Ez = Z ■ 4.5 PeV. The cut-off for the 
heaviest elements agrees with the energy of the second knee at ~ 400 PeV, which is 
interpreted as the end of the galactic CRs, while the knee is caused by the cut-off 
of the light elements. The sum spectrum of all elements is given by the solid line, 
which fits nicely the average measured spectrum up to 100 PeV. At low energies 
where the nuclei traverse a large amount of matter (~ 10 g/cm 2 ), heavy nuclei are 
more likely to interact with the interstellar matter as compared to light elements 
(iiinei oc A 2 / 3 ) and the spectra observed at earth are expected to be slightly flatter 
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Fig. 7. Mean logarithmic mass as function of energy obtained by direct observations (dark grey 
area) and air shower exDeriments (light grey area) coniDarcd with different models (line s') . a) Ac- 
celeration in 

SNR J2()|22|21|2;H 

b) acceleration in GRBs 4y ^ 4 *^ t> , single source model2^ reaccclcr- 
ation in the galactic win c) diffusion in G alay l ^ l ^ l ^ l d - ) pro pagation in theGalaxjEEH 
as well as interaction with background photonJ^Sl and neutrinoJSSJ For details se d^7 | 



for heavy nuclei. At the respective knees A is less than 1 g/cm 2 , thus for the 
heaviest elements around 400 PeV more than 40% of the nuclei are expected to 
survive without interactiorP^. 



5. THE KNEE IN THE ENERGY SPECTRUM 

The bulk of CRs is assumed to be accelerated in strong shock fronts of SNrJ^HI. The 
finite lifetime of a shock front (~ 10 5 a) limits the maximum energy attainable for 
particles with charg e Z to Em „, x ~ Z ■ (0.1 — 5) PeV. Many versions of this scenario 
have been discussecP^^^222T The models differ in assumptions of properties of 
the SNRs like magnetic field strength, available energy etc. This yields differences 
in (In A), as can be inferred from Fig.[7]i. While older mo dela^H limit the maximum 
energy to about 0.1 PeV, recent idea*^ 1 , taking into account latest observations of 
SNRs, predict maximum energies above 1 PeV. In such a model sufficient energy 
is released from SNRs to explain the observed spectra, see in Fig. [21 the calcula- 
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tions by Sveshnikova et ai. A special case ol SNR acceleration is the single source 
model , which predicts in the knee region pronounced structures in the all-particle 
energy spectrum, caused by a single SNR. Such structures can not be seen in the 
compilation ol Fig. El 

In the literature also other accele ration m echanisms , like the acceleration ol 
particles in 7-ray bursts, are They differ in their interpretation 

ol the origin lor the knee. The approach by Plaga, assuming Fermi acceleration in 
a "cannon ball" is not compatible with the measured (In A) values, see Fig. 03. 
A different interpretation ol acceleration in the cannon ball model yields - at the 
source - a cut-off lor individual elements proportional to their mass due to effects 
ol relativistic beaming in jets. The predictions ol the actual model are compatible 
with recent datcP^. However, it remains to be clarified how a detailed consideration 
ol the propagation processes, e.g., in a diffusion model, effects the cut-off behavior 
observed at earth. Gamma-ray bursts as a special case ol supernova explosions are 
proposed^! to accelerate CRs from 0.1 PeV up to the highest energies (> 10 20 eV). 
In this approach the propagation ol CRs is taken into account and the knee is 
caused by leakage from the Galaxy leading to a rigidity dependent cut-off behavior. 

The propagation is accompanied by leakage ol particles from the Galaxy. With 
increasing energy it becomes more and more difficult to confine the nuclei to the 
Galaxy. As mentioned above, the pathlength decreases as A oc E~ 5 . Such a decrease 
will ultimately lead to a complete loss of the particles, with a rigidity dependent 
cut-off of the flux for individual ele ments. Many ap proaches have been un dert aken 
to describe the propagation procesJ ^^^^^ I The Leaky Box modefSl and 
the anomalous diffusion modei^ yield cut-offs significantly weaker than the data 
shown in Fig. [21 

The propagation as described in diffusion modelJ ^^^ I yields (lnyl)-values 
which are presented in Fig.U|;. The models are based on the same principal idecPl 
but take into account different assumptions on details ol the propagation process, 
like the structure ol galactic magnetic fields etc. This results in a more or less strong 
cut-off for the flux at the individual knees and, ac cor dingly, in a more or less strong 
increase of (In A) . The model by Kalmykov et alW^ has been used to describe the 
observed spectra in Fig. 

During the propagation pha se, reacceleration ol particles has been suggested at 
shock fronts in the galactic wino™4 Also this mechanism yields a rigidity dependent 
cut-off. 

Another hypothetical explanation for the knee are interactions ol CRs with back- 
grou nd par ticles like massive neutrino J5EEH1 or photo disintegration in dense photon 
field er" ' ' ! ^ ! Such models appear to be excluded with a high level ol confidence. The 
interactions would produce a large amount ol secondary protons, which results in 
a light mass composition at high energies, not observ ed by t he experiments, see 
Fig. 0i. Furthermore, a massive neutrino, proposed in I I can be excluded by 
measurements ol the WMAP and 2dFGRS experiments^. 

A completely different reason for the knee is the idea to transler energy in 
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nucle on-n ucleon interactions into particles, like gravitonJ^or extremely high-energy 
muonJ^, which arc not observable (or not yet observed) in air shower experiments. 
The latter pro posal seems to be excluded by recent measurements of the Baikal 
experiment!^ setting upper limits for the flux of muons above 10 5 GeV. 

6. CONCLUSION 

During the last decade significant progress has been made in the measurement of 
galactic CRs. Summarizing the large number of experimental observations, there 
are indications for a standard picture. At least a large fraction of CRs seems to 
be accelerated in supernova remnants up to energies of Z ■ (0.1 — 5) PeV. Higher 
energies may be reached in additional sources, such as 7-ray bursts. The elemental 
composition of the accelerated material is extremely similar to that in the solar 
system. The particles propagate in a diffusive process through the Galaxy. With 
rising energy the pathlength decreases and the particles escape easier from the 
Galaxy. This brings about the knee in the energy spectrum. The general shape of 
the energy spectra should be determined by the propagation process, maybe slightly 
modulated by properties of the source spectra. 
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